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In this article mixed CuInP2(S I Sei_ ;E )6 crystals were investigated by broadband dielectric spec- 
troscopy (20 Hz - 3 GHz). The complete phase diagram has been obtained from these results. The 
phase diagram of investigated crystals is strongly asymmetric - the decreasing of ferroelectric phase 
transition temperatures in CuInP2(S a; Sei_ : r)6 is much more flat with small admixture of sulphur 
then with small admixture of selenium. In the middle part of the phase diagram (x=0.4-0.9) the 
dipolar glass phase has been observed. In boundary region between ferroelectric order and dipolar 
glass disorder with small amount of sulphur (x=0.2-0.25) at low temperatures the nonergodic relaxor 
phase appears. The phase diagram was discussed in terms of random bonds and random fields. 

PACS numbers: 77.22.-d, 77.80.-c, 77.22. Gm, 81.30. -t 



I. INTRODUCTION 

Solid systems present many interesting types of phase 
transitions, with ferro, antiferro, or modulated long range 
order at lower temperatures. Disordered cooperative sys- 
tems have also attracted a lot of attention. Nonergodic 
relaxor, dipolar glass phases or coexistence of ferroelec- 
tric and dipolar glass phases can appear in disordered 
systems at low temperatures. The nature of these phases 
continues to generate considerably experimental and the- 
oretical interest. 

CuInP2S6 crystals represent an unusual example of 
a anticollinear two-sublattice ferriclcctric system [1, 2, 
3, 4]. Here a first-order phase transition of the order- 
disorder type from the paraelectric to the ferriclcctric 
phase is realized (T c = 315 K). The symmetry reduction 
at the phase transition (C2/c — > Cc) occurs due to order- 
ing in the copper sublattice and displacement of cations 
from the centrosymmetric positions in the indium sublat- 
tice. The spontaneous polarization arising at the phase 
transition to the ferrielectric phase is perpendicular to 
the layer planes. These thiophosphatcs consist of lamel- 
lae defined by a sulphur framework in which the metal 
cations and P - P pairs fill the octahedral voids; within 
a layer, the Cu, In, and P-P form triangular patterns 
[1, 2, 3]. The cation off-centering, 1.6 A for Cu 1 and 0.2 
A for In 111 , may be attributed to a second-order Jahn- 
Teller instability associated with the d 10 electronic con- 
figuration. The lamellar matrix absorbs the structural 
deformations via the flexible P2S6 groups while restrict- 
ing the cations to antiparallcl displacements that mini- 
mize the energy costs of dipole ordering. Each Cu ion 
can occupy two different positions. The Cu, In and P-P 
form triangular patterns within the layer. Rclaxational 
rather than resonant behaviour is indicated by the tem- 
perature dependence of the spectral characteristics, is in 
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agreement with X— ray investigations. It was suggested 
that a coupling between P2S 6 deformation modes and 
Cu 1 vibrations enables the copper ion hopping motions 
that lead to the loss of polarity and the onset of ionic 
conductivity in this material at higher temperatures [4]. 
The investigation of ionic conductivity in CUI11P2S6 [5, 6] 
have showed that croc follows the Arrhcnius law with the 
activation energy Ea = 0.73 cV [5] and more detailed in- 
vestigations showed Ea=0-&55 eV [6]. 

The results of dielectric investigations of Culid^Scg 
showed two phase transition: a second-order one at 
Ti= 248 K and a first-order transition at T c =236 K 
[ ]. The results followed to the conclusion that an in- 
commensurate phase occurs between Ti and T c . How- 
ever, the calorimctric investigations showed only a broad 
phase transition between 220 and 240 K in this com- 
pound [8] . More accurate broadband dielectric investiga- 
tions showed only nearly second order phase transition 
at T" c =226 K [9]. From a single-crystal X-ray diffrac- 
tion study follows that the high- and low-temperature 
structures of CuInP2Sc6 (trigonal space group P-31c and 
P31c, respectively) are very similar to those of CuInP 2 S6 
in the paraelectric and ferrielectric phases, with the Cu 7 
off-centering shift being smaller in the former than in the 
latter [3, 8]. There the thermal evolution of the cell pa- 
rameters of CuInP 2 Se 6 was obtained by full profile fits 
to the X-ray diffractograms. Both cell parameters a and 
c slightly decrease on cooling, and a parameter shows 
a local minimum at T=226 K. This behaviour is quite 
different from the anomalous increases found in the cell 
parameters of CuInP2Sg when heating through the tran- 
sition [1, 3]. 

The important feature of sclcnidcs is the higher co- 
valcncc degree of their bonds. Evidently, for this rea- 
son the copper ion sites in the low-temperature phase 
of CuInP2Se6 are displaced only by 1.17 A [8] from the 
middle of the structure layers in comparison with the cor- 
responding displacement 1.6 A for CuInP2S6 [1]. These 
facts enable to assume that the potential relief for cop- 
per ions in CuInP2Se6 is shallower than for its sulphide 
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analog. Presumably, for this reason the structural phase 
transition in the selenide compound is observed at lower 
temperature than for the sulphide compound. Prelimi- 
nary dielectric investigations of CuInP2(S ;c Sei_ :1 -)6 crys- 
tals are presented in [11, 12]. The data in [11] is measured 
only at frequency 10 kHz and the paper [12] contain data 
only on one compound - CuInP2 (So. 7800.3)6- 

The aim of this paper is to investigate phase diagram 
of mixed CuInP2(S 2; Sei_ :E )6 crystals via broadband di- 
electric spectroscopy. We showed that in mixed crystals 
with the increasing amount of impurities two smearing of 
ferroelectric phase transition scenarios are possible: fer- 
roelectric - inhomogencous ferroelectric - dipolar glass or 
ferroelectric - relaxor - dipolar glass. 



II. EXPERIMENTAL 

Crystals of CuInP2(S ;c Sei_ :1 -)6 were grown by Bridg- 
man method. For the dielectric spectroscopy the plate 
like crystals were used. All measurements were per- 
formed in direction perpendicular to the layers. The 
complex dielectric permittivity e* was measured using 
the HP4284A capacitance bridge in the frequency range 
20 Hz to 1 MHz. In the frequency region from 1 MHz to 3 
GHz measurements were performed by a coaxial dielec- 
tric spectrometer with vector network analyzer Agilent 
8714ET. All measurements have been performed on cool- 
ing with controlled temperature rate 0.25 K/min. Silver 
paste has been used for contacting. 



III. 



RESULTS AND DISCUSSION 



A. Influence of small amount of sulphur to phase 
transition dynamics in CuInP2See crystals 

A small amount of admixture can significant changes 
properties of ferroelectrics. In mixed CuInP2(S 2; Sei_ a ;)6 
crystals with a;<0.1 the ferroelectric phase transition is 
observed (Fig. 1). Here the dielectric permittivity max- 
imum temperature (T m ) is frequency-dependent only at 
higher frequencies (above 1 MHz). The phase transition 
temperature can be defined by T m at low frequencies (be- 
low 1 MHz) . The temperature behaviour of the dielectric 
dispersion of CUI11P2SC6 crystals with a small admixture 
of sulphur (Fig. 2) is very similar to the dielectric disper- 
sion of pure CuInP2Seg crystals [9]. At higher tempera- 
tures (T >> T c ) the dielectric dispersion reveals in 10 8 
- 10 10 Hz frequency range. With decreasing temperature 
the dielectric dispersion become broader and appears at 
lower frequencies. At lower temperatures (T << T c ) the 
dielectric dispersion remains in the 10 6 - 10 10 Hz fre- 
quency range and only its strength decreases on cooling. 

More information about the phase transition dynamics 
can be obtained by analysis the dielectric dispersion with 



the Cole- Cole formula 

e*(y) = Soo + 



Ae 



1 + {iujTccY~ acc 



(1) 



where Ae represents dielectric strength of the relaxation, 
tcc is the mean Cole-Cole relaxation time, £oo repre- 
sents the contribution of all polar phonons and elec- 
tronic polarization to the dielectric permittivity and acc 
is the Cole-Cole relaxation time distribution parameter; 
when acc=®i Eq. 1 reduces to the Dcbye formula. Ob- 
tained parameters are presented in Fig. 3. The Cole- 
Cole parameters of all presented compounds show the 
similar behaviour: the Cole-Cole distribution parame- 
ter acc strongly increases on cooling, reciprocal dielec- 
tric strength 1 / Ae exhibits a minimum at ferroelectric 
phase transition temperature, the soft mode frequency 
v r = 1/(2tttcc) slows down 011 cooling in the paraelec- 
tric phase. The temperature dependence of the dielectric 
strength Ae was fitted with the Curie- Weiss law (Fig. 3) 



Ae = C p j/(\T-T c \), 



(2) 



where C p j is the Curie-Weiss constant and Tc is the 
Curie- Weiss temperature. The temperature dependence 
of soft mode frequency v r in paraelectric phase was fitted 
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FIG. 1: Temperature dependence of the complex dielectric 
permittivity of CuInP2(So.iSeo.g)6 crystals measured at sev- 
eral frequencies. 
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with the equation 



v r =A{T-T c ), 



(3) 



where A is a constant. Obtained parameters are pre- 
sented in Table 1. The phase transition temperature Tq 
in mixed crystals strongly decreases from 225 K to 185 
K. For all the compounds the C p /Cf ratio is about 1.5, 
for the second order phase transitions this ratio must be 
2, for the first order one - higher than 2. The assumption 
was made that in these crystals between paraelectric and 
ferroelectric phase an additional incommensurate phase 
exists [11]. However, in all mixed CuInP 2 (S K Sei_ :i; ) 6 
crystals with x<0.1 no anomaly above the main (ferro- 
electric) phase transition was observed (Fig. 1). 

Below the ferroelectric phase transition temperature 
the dielectric dispersion is broad and part of it appears 
in the low frequency region (Fig. 1). This part is caused 
by ferroelectric domain dynamics. Therefore, the contri- 
bution of ferroelectric domain dynamics effectively raises 
the dielectric strength Ae in the ferroelectric phase and 
Cf constant. 
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FIG. 2: Frequency dependence of the complex dielectric per- 
mittivity of CuInP2(So.iSeo.9)6 crystals measured at several 
temperatures. Lines are results of Cole-Cole fits. 



TABLE I: Parameters of phase transition dynamic of 
CuInP 2 Seg crystals with small admixture of sulphur (s<0.1). 



compound 


C p , K 


C p / C f 


A, MHz/K 


T c , K 


CuInP 2 Se 6 from [!)] 


591.7 


1.33 


271.9 


225 


CuInP 2 (Se .98So.o2)6 


309.6 


1.43 


193.4 


215.7 


CuInP 2 (Se .95So.o5)6 


980.3 


1.66 


79.3 


208.2 


CuInP 2 (Se . 9 So.i)6 


2380.9 


1.52 


44.4 


185 
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FIG. 3: Temperature dependence of the Cole-Cole 
parameters of complex dielectric permittivity for the 
CuInP 2 (S a; Sei_ a ;)6 crystals with a;<0.1. The v r lines were 
obtained from fit with Eq. 3 and the 1/Ae lines were ob- 
tained from Curie- Weiss fit. The data for CuInP 2 Se6 is from 
[9]. 



B. Nonergodic relaxor phase in mixed 
CuInP 2 (S a; Sei_ :E )6 crystals 

Recently, the relaxor-like behaviour as an embryo of 
the glass state is proposed in the antiferroelectric-glass 
phase boundary region of DRADP crystals family [18]. 
Here it is showed that the growth of glass ordering is 
in quite a different pattern from that of the ferroelectric- 
glass phase boundary region. In this section we shall pre- 
sented two very similar CuInP2(S 2; Sei_ 2: )6 compounds 
(x=0.2 and x=0.25), which exhibit peculiar dielectric be- 
haviour. Each composition shows just one maximum in 
e'(T) and e"(T) in the range of 110 and 145 K at frc- 
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quency 10 kHz [11]. The temperature dependences of the 
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FIG. 4: Temperature dependence of the complex dielectric 
permittivity of CuInP2(So.25Seo.7s)6 crystals measured at sev- 
eral frequencies. 

complex dielectric permittivity e* at various frequencies 
of these crystals show typical relaxor behaviour. As an 
example, dielectric permittivity of CuInP2(Seo. 7580.25)6 
crystal is shown in Fig. 4. There is a broad peak in the 
real part of dielectric permittivity is observed. With fre- 
quency T m and the magnitude of the peak increases in 
the whole frequency range. There is a strong dielectric 
dispersion in a radio frequency region around and below 
T m at 1 kHz. The value of T mm (the temperature of 
the maximum of losses) is much lower than that of T m 
at the same frequency. The position of the maximum of 
dielectric permittivity is strongly frequency-dependent; 
no certain static dielectric permittivity can be obtained 
below and around dielectric permittivity maximum tem- 
perature T m at 1 kHz. Such behaviour can be described 
by the Vogel-Fulcher relationship 

Ef 



vq exp 



k(T m - T ot ) ' 



(4) 



where k is the Boltzman constant, Ef, i>q, T ot are pa- 
rameters of this equations. Obtained parameters are pre- 
sented in Table II. 
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FIG. 5: Frequency dependence of the complex dielectric per- 
mittivity of CuInP2(So.25Seo.75)6 crystals at several tempera- 
tures. Lines are results of fits with distributions of relaxation 
times (solid) and of Cole-Cole fit (dot). 



TABLE II: Parameters of the Vogel-Fulcher fit of the T m de- 
pendence of frequency for CuInP2(S I Sei_ a; )6 crystals with 
0.2<z<0.25. 



compound 


uo, GHz 


Tot, K 


E f /k, K 


CuInP 2 (Seo.75So.25)6 


38.34 


96.8 


370 


CuInP2(Seo.8So.2)6 


10.96 


134.5 


150 



The dielectric dispersion of CuInP2(Seo. 7580.25)6 crys- 
tals show strong temperature dependence (Fig. 5): at 
higher temperatures the dielectric dispersion is only in 
10' - 10 10 Hz region, on cooling the dielectric dispersion 
becomes broader and more asymmetric. Strongly asym- 
metric and very broad dielectric dispersion is observed 
below dielectric permittivity maximum temperature T m 
at 1 kHz. The Cole-Cole formula (Eq. 1) can describe 
such dielectric dispersion only at higher temperatures, 
due to predefined symmetric shape of the distribution of 
the relaxations times. This is clearly visible in Fig. 5, 
where the Cole-Cole fit is shown as doted line. Not only 
Cole-Cole formula, however, other very well known pre- 
defined dielectric dispersion formulas, such as Havriliak- 
Negami, Cole-Davidson cannot adequate describe the di- 
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FIG. 6: Relaxation time distribution 

CuInP2(So.25Seo.75)6 crystals at various temperatures. 
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electric dispersion of the presented crystals. More general 
approach must be used for determination of the broad 
continuous distribution function of relaxation times /(r) 
by solving a Fredholm integral equations 



s'(lo) 



As 



f(r)d(lnT) 



Loo 1+W 2 T 



1 + L0 2 T 2 

LOT f (T)d(lriT) 
2 



with the normalization condition 



f(T)d(ln T ) = 1. 



(5a) 
(5b) 

(0) 



The most general method for the solution is the 
Tikhonov regularization [14, 15] method. The calculated 
distribution of relaxation times of CuInP 2 (So. 2 5Seo.75)6 
crystals is presented in Fig. 6. The symmetric and 
narrow distribution is observed only at higher tempera- 
ture T >> T m (at 1 kHz), on cooling the distributions 
becomes broader and more asymmetric so that be- 
low T m (at 1 kHz) second maximum appears. Such 
behaviour of distribution of relaxation times have 
been already observed in a very well known relaxors: 
Pb(Mg 1/3 Nb 2/3 )03 (PMN) [10], Pb(Mg 1/3 Nb 2/3 )0 3 - 
Pb(Zn 1/3 Nb 2/3 )0 3 -Pb(Sc 1/2 Nb 1/2 )0 3 (PMN-PZN- 
PSN) [17], Pb(Mg 1/3 Ta 2/3 )0 3 (PMT) [18] and 
Sr . 6 iBa . 39 Nb 2 O 6 (SBN) [19]. From calculated 
distributions of relaxation times the most probable 
relaxation time r mp , longest relaxation time T max 
and T rnin shortest relaxation time (the level 0.1 was 
chosen as sufficient accurate) has been obtained (Fig. 



TABLE III: Parameters of the Vogel-Fulcher fit of the tem- 
perature dependencies of the longest relaxation times T ma x m 
CuInP 2 (S a: Sei- 2 ,)6 crystals with 0.2<x<0.25. 



compound 




To, K 




CuInP 2 (Se .75So.25)6 


2.52xl(T 8 


118.9 


60.5 


CuInP 2 (Se . 8 So.2)6 


1.02X10 -10 


129.4 


211.01 



7). The shortest relaxation time T m i n is about 0.1 
ns for CuInP 2 (So. 2 5Se .75)6 and about 0.01 ns for 
CuInP 2 (So. 2 Seo.8)6; it increases slowly with the increase 
of temperature. The longest relaxation time r max 
diverges according to the Vogel-Fulcher law 



T~0r, 



; exp 



k(T-T )' 



(7) 



where To is the freezing temperature, E max is the acti- 
vation energy of the longest relaxation times r max and 
T 0max is the longest relaxation time at very high tem- 
peratures. The obtained parameters are presented in Ta- 
ble III, however the most probable relaxation time r mp 
diverges with good accuracy according to the Arrhcnius 
law: 



E 

^mp 

T~mp = TOmp exp ^ . 



(8) 



where E mp is the activation energy of the most proba- 
ble relaxation times T mp , and TQ mp is the most probable 
relaxation time at very high temperatures. Obtained pa- 
rameters are To„ lp = 4.6xl0~ 16 s and i? mp /k = 2365.3 
K for CuInP 2 (Seo.75So. 2 5)6 and r 0mp = 1.2xl0~ 14 s and 
E mp /k = 1806.3 K for CuInP 2 (Seo. 8 So. 2 )6- Such phe- 
nomenon can be caused by a distribution of Vogel-Fulcher 
temperatures T , where 0<r o <T maa:; o [20], [21]. In our 
case T max would correspond to a Vogel-Fulcher tempera- 
ture of T max and is the freezing temperature of the most 
probable relaxation time and all shorter relaxation times. 
The temperature dependence of the reciprocal static di- 
electric permittivity l/e(0) was fitted with sperical ran- 
dom bond random field (SRBRF) 



6(0) = 



C p (l - q E A) 
kT - J(l - q E A) 



(9) 



where J is the mean coupling constant and q E A is 
Edwards- Anderson order parameter, if q E A=0 then this 
equation becomes the Curie- Weiss law. The Edwards- 
Anderson order parameter q E A for relaxor can be deter- 
mined by equation [22]: 

qEA = (^) 2 (q E A + J^)(l-qEA) 2 , (10) 

where A J is the variance of the coupling and A/ is 
the variance of the random fields. Obtained parameters 
we will discussed further together with random bonds 
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FIG. 7: Temperature dependence of the longest r max , most 
probable r mp , shortest T m i n relaxation times and static di- 
electric permittivity e(0) in CuInP2(S a; Sei_ a ;)6 crystals, with 
x=0.2 (solid points) and £=0.25 (open points). The t(T) 
lines were obtained from Vogel-Fulcher (for longest relaxation 
times) and from Arrenius (for most probable relaxation times) 
fits. The static dielectric permittivity e(0) lines were obtained 
from Eqs. 9 and 10. 

random fields parameters of other mixed crystals. We 
must admit that the equations of the SRBRF model de- 
scribe well static dielectric properties of the presented 
crystals. At sulphur concentrations between £=0.25 and 
2=0.2, morphotropic phase boundary between the para- 
electric phases C2/c (characteristic for CuInP 2 S6) and 
P-31c (characteristic for CuIuE^Seg) or respectively fer- 
rielectric phases Cc and P31c were suggested [11]. These 
results were later confirmed by X-ray and Raman investi- 
gations [23]. Therefore, the disorder in these mixed crys- 
tals is very high, and it can be reason of relaxor nature 
of the presented crystals. 

C. Dipolar glass phase in mixed CuInP2(S a ,Sei_ :c )6 
crystals 

For CuInP2(S- E Sei_ x )6 crystals with x=0.4-0.9 no 
anomaly in static dielectric permittivity indicating the 
polar phase transition can be detected down to the lowest 
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FIG. 8: Temperature dependence of the complex dielectric 
permittivity of CuInP2(So.sSeo.2)6 crystals measured at sev- 
eral frequencies. 



temperatures. The dielectric spectra of these crystals are 
very similar. As an example, real and imaginary parts of 
the complex dielectric permittivity of CuInP2(So.sSeo.2)6 
crystals are shown in Fig. 8 as a function of temperature 
at several frequencies. It is easy to see a broad disper- 
sion of the complex dielectric permittivity starting from 
260 K and extending to the lowest temperatures. The 
maximum of the real part of dielectric permittivity shifts 
to higher temperatures with increase of the frequency to- 
gether with the maximum of the imaginary part and man- 
ifests typical behaviour of dipolar glasses. The dielectric 
dispersion is symmetric of all crystals under study so that 
it can easily be described by the Cole-Cole formula (Fig. 
9). The temperature dependence of the Cole-Cole pa- 
rameters confirms typical behaviour for dipolar glasses 
(Fig. 10): the mean Cole-Cole relaxation time diverge ac- 
cording to the Vogel-Fulcher law (Eq. 7), the Cole-Cole 
distribution parameter olqc strongly increases on cooling 
and reaches value 0.5 below 100 K, the static dielectric 
permittivity temperature dependence has no expressed 
maxima. Usually such behaviour is analyzed in terms 
of the three-dimensional random-bond random-field (3D 
RBRF) Ising model of Pirc et al [24]. In terms of this 
model, the temperature dependence of static dielectric 
permittivity can be described with the Eq. 9. The order 
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FIG. 9: Frequency dependence of the complex dielectric per- 
mittivity of CuInP2(So.8Seo.2)6 crystals at several tempera- 
tures. Lines are results of Cole-Cole fits. 



parameter is defined by the two coupled self-consistent 
equations[25] 

/°° dz , , 77 . , z 2 . , 
- — —tanh( — )exp( ), (11 



QEA = J ^ J^oJ tanh2 (kf) ex P(~Y^ (12) 
where P is the polarization and 

7] = (AJ 2 q EA + Aff 5 Z + JP. (13) 

The Equation 9 describe good enough static dielec- 
tric properties of presented dipolar glasses and obtained 
parameters are in good agreement with parameters ob- 
tained from Vogel-Fulcher fits, according to formula [26] 

T a = AJ/k B . (14) 

Obtained parameters we will discuss further below to- 
gether with random bonds random fields parameters of 
other mixed crystals. 



FIG. 10: Temperature dependence of the Cole-Cole 
parameters of complex dielectric permittivity for the 
CuInP2(Sa;Sei-a;)6 crystals with 0.4<x<0.9. The r lines were 
obtained from Vogel-Fulcher fit and the e(0) lines were ob- 
tained from 3D RBRF model fit. 



D. Influence of small amount of selenium to phase 
transition dynamics in CuInP2Ss crystals 

Temperature dependence of the dielectric permittiv- 
ity of CuInP2S6 crystals with a small amount of sele- 
nium (x=0.98) is presented in Fig. 11. A small amount 
of selenium changes dielectric properties of CuInP 2 S6 
crystals significantly: the temperature of the main di- 
electric anomaly shift from about 315 to 289 K, the 
maximum value of the dielectric permittivity e' signifi- 
cantly decreases from about 180 to 40 (at 1 MHz), at 
higher frequencies (from about 10 MHz) the peak of 
dielectric permittivity becomes frequency- dependent in 
CuInP2(So.98Serj.rj2)6 crystals and a critical slowing down 
disappears [6]. An additional dielectric dispersion ap- 
pears at low frequencies and at low temperatures. The 
CuInP2(So.g5Seo.o5)6 crystals exhibit qualitatively simi- 
lar dielectric anomaly with T c and e' max shifting to lower 
values. The dielectric dispersion of presented crystals is 
symmetric (Fig. 12) so that it can be correctly described 
by the Cole-Cole formula (Eq. 1). The Cole- Cole param- 
eters are shown in Fig. 13. The parameters of the Cole- 
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FIG. 11: Temperature dependence of the complex dielectric 
permittivity of CuInP2(So.9sSeo.o2)6 crystals measured at sev- 
eral frequencies. 



TABLE IV: Parameters of phase transition dynamic of 
CuInP2S6 crystals with small admixture of selenium. 



compound 


C P ,K 


C p /Cf 


T Cp , K 


Tc/,K 


CuInP2(Seo.o5So.95)6 


8587.7 


2.99 


137.2 


368.7 


CuInP 2 (Seo.o2S .98)6 


1906.5 


7.01 


236.9 


282.6 
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FIG. 12: Frequency dependence of the complex dielectric per- 
mittivity of CuInP2(So.g8Seo.o2)6 crystals measured at several 
temperatures. Lines are results of Cole-Cole fits. 



TABLE V: Parameters of the Vogel-Fulcher fit of the tem- 
perature dependencies of the mean relaxation times tcc in 
CuInP2(S I Sei_ I )6 inhomogeneous ferroelectrics. 

compound tq, s To, K E/k, K 



Cole distribution of relaxation acc strongly increase on 
cooling and reach 0.43 at low temperatures. The temper- 
ature dependence of the dielectric strength Ae was fitted 
with the Curie- Weiss law (Eq. 2). Obtained parameters 
are summarized in Table IV. The difference Tc P -Tcf 
and ratio C p /C/ in these crystals indicate a first order, 
order-disorder phase transition. In ferroelectric phase the 
mean relaxation time Tcc decreases only in a narrow 
temperature region and only for CuInP2(So.9sSeo. 02)61 
further on cooling a significant increasing of times tcc is 
observed. This increasing can be easily explained by the 
Fogcl-Vulchcr law (Eq. 7). These parameters arc sum- 
marized in Table V. Note that all parameters of different 
compounds in Table V are close to each other. Such a 
behaviour is very similar to behaviour of bctainc phos- 
phite with a small amount of betaine phosphate [27] and 
in RADA [28] crystals, where a proposition that a coexis- 



CuInP2(Se .9 5 So.o5)6 8.5xHT 12 1150 31 
CuInP 2 (Seo.98So.o2)6 3.77X10" 11 1215 28 



tence of the ferroelectric order and dipolar glass disorder 
appears at low temperatures was proposed. Therefore we 
can conclude that mixed Cubd^S^Sei-^g crystals with 
x>0.95 also exhibit at low temperatures a coexistence of 
ferroelectric and dipolar glass disorder. 



E. Phase diagram 

In this section we will discuss phase diagram in terms 
of random bonds and random fields. For ferroelectrics we 
assume that mean coupling constant J/k is equal to T@, 
because Curie- Weiss fit is accurate for these compounds 
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Temperature (K) Temperature (K) Temperature (K) 



FIG. 13: Temperature dependence of the Cole-Cole 
parameters of complex dielectric permittivity for the 
CuInP2(S a; Sei_ a ;)6 crystals with x=0.95 (open points) and 
x=0.98 (solid points). The r lines were obtained from Vogel- 
Fulcher fit and the Ae lines were obtained from Curie- Weiss 
fit. 



and in this case Eq. 9 becomes Curie- Weiss law. Also for 
crystals with with x<0.1, for the same reason we assume 
that AJ and A/ are 0. For ferroelectrics with x>0.95 we 
obtained AJ from T (Eq. 14), we assumed that A/=0. 

In Fig. 14 we present the obtained phase diagram of 
mixed crystals. In the mixed CuInP2(S 2; Sci_ 2 ;)6 with 
x>0.95 and x<0.1 crystals the mean coupling constant 
J>(A/+A J 2 ) 05 , therefore, they undergo ferroelectric 
phase transition at J /k. However is significant difference 
between phase transition dynamics of mixed crystals with 
s>0.95 and a;<0.1. In mixed crystals with a;<0.1 no any 
coexistence of ferroelectric order and dipolar glass disor- 
der is observed down to the lowest temperature (80 K). 
At temperatures below 100 K the dielectric permittivity 
of these compounds is very low (about 3), therefore, the 
phase coexistence in these compounds is unlikely. In the 
ferroelectric phase these crystals split into domains, it is 
evidenced in low frequency dielectric dispersion spectra 
(Fig.l). However similar ferroelectric domains already 
are observed in pure CuInP 2 Se 6 crystals [9]. Really, in- 
fluence of small amount of sulphur to phase transition 
dynamics of mixed crystals appears only by reduction 




0.0 0.2 0.4 0.6 0.8 1.0 



sulphur concentration 



FIG. 14: Phase diagram of the mixed CuInP2(S I Sei_ a: )6 
crystals (FE - ferroelectric phase, NR - nonergodic relaxor 
phase, DG - dipolar glass phase, FE+DG - ferroelectric and 
dipolar glass coexistence). 

Tc (Table I). The influence of small amount of sele- 
nium to phase transition dynamics is more significant - 
already at x=0.95 the ferroelectric phase transition in 
Tec is less expressed (Fig. 13). Such influence is ex- 
pressed also in other properties: rapid decreasing in Tc, 
appearance of ferroelectric and dipolar glass phase coex- 
istence at £=0.98 and onset of dipolar glass disorder with 
x between 0.9 and 0.95. 

For crystals with x=0.2 and 0.25 J<(A/+A J 2 ) 5 and 
J~ (A/+A J 2 ) 5 therefore the nonergodic relaxor phase 
appears in these crystals at low temperatures. In the 
presence of an external electric field E meaning coupling 
constant J is expected to vary as 

J(E) = J(0) + aE 2 . (15) 

For electrical field E that J{E) > (A/+A J 2 ) 5 , in mixed 
crystals should be observed relaxor to ferroelectric phase 
transition. The possible existence of relaxor phase in 
mixed ferroelectric-antiferroelectric crystals is stated in 
[13, 29]. Really, no any evidence is indicated for po- 
lar nanoregions existence in mixed crystals. We try to 
fill this gap of information presenting two mixed crys- 
tals, where dielectric behaviour is very similar to very 
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well known relaxors PMN [30] and SBN [31] (the differ- 
encies are only in T rn and e' max values). On the other 
hand, in phase diagram with less selenium concentration 
no area with nonergodic relaxor phase (Fig. 14) appears. 
The main cause of such phase diagram is that disorder 
((A/+AJ 2 ) - 5 ) is highest at 2=0.2, where mean coupling 
constant is also high enough. Usually, for mixed crystals 
is assumed that concentration dependence for A/ is such 
[24] 

Af = 4x(l - x)Af max . (16) 

For AJ similar behaviour also was assumed. In this case 
if J has minimum at x=0.5 the nonergodic relaxor phase 
can not be observed. However any existing theories can 
not explain A J and A/ concentration dependence. 

For compounds 0.9>x>0.4 the relation J<s; 
(A/+AJ 2 ) 5 is valid, consequently in these compounds 
a dipolar glass phase appears at low temperatures. 



By further increasing selenium concentration the dipolar 
glass phase appears. In contrast to in CuIut^Seg even a 
high concentracion of admixture of sulphur (x=0.1) has 
no any influence to the feroelectric order. The some 
degree of ferroelectric order exist even for 1=0.2 and 
2=0.25, however, in these crystals the ferroelectricity is 
broken into polar nano regions. The random bonds and 
random fields model clearly describe the asymmctricity of 
phase diagram of mixed CuInP2(S x Sei_ x )6, however this 
model can not identified origin of the effect. To summa- 
rize, the first experimental evidence for smearing noner- 
godic relaxor phase into dipolar glass phase by some dop- 
ing is presented. For other relaxors the search of some ad- 
mixture which transforms relaxor state into dipolar glass 
can also be performed. 



IV. CONCLUSIONS 



The ferroelectric order in CUI11P2S6 is reduced already 
for small (x=0.98) substitution of sulphur by selenium. 
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